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TECHNICAL PAPER 


EFFECT OF TYPE OF LOAD ON STRESS ANALYSIS 
OF THIN- WALLED DUCTS 

I. INTRODUCTION 


The standard procedure for qualifying the design of duct (pipe) systems in the space shuttle 
main engine (SSME) has been fairly well defined. However, since pipe elbows are quite common and 
important in the SSME duct systems, a clear understanding of the detailed stress profile of these 
components is necessary for accurate structural and life assessments. The analytical work on the 
analysis of pipe bends was originated by Von Karman 1 and continued by several researchers 2 3 for 
decades. In recent years, notable contributions on this topic have been made by Spence and 
Thompson, 4 Thomson and Spence, 5 Whatham and Thompson, 6 Whatham, 7 Karbin, Rodabaugh, and 
Whatham, 8 and Natarajan, 9 using the shell theory and finite element method (FEM). However, all 
these analyses provide limited information about the stress distributions at the interface of the 
elbow with the straight pipes (tangent point). In general, the weld location at the tangent point is the 
most critical location governing the design of duct systems. This is especially true in the design of 
SSME duct systems because of the harsh environment (high pressure/temperature), large dynamic 
loads, and weld discontinuity (mismatch/weld bead). A space vehicle design also requires 
optimization of weight for maximum lift-off capability, thus resulting in thin-walled duct systems. 
Additionally, hot fire data from the strain gauges mounted in the vicinity of these welds are used to 
anchor these data to the dynamic models for predicting the loads at other locations. Based on the 
previous studies mentioned above, it is well known that the circumferential and longitudinal stresses 
at the cross section of the elbow midsection do not peak at the same location. Furthermore, the 
location (position angles of cross section, fig. 1) of peak stress of these stresses is not at the 
expected position angle. Therefore, this study was initiated to predict the stress profile at/near the 
tangent point along the cross section of the duct under various types of loads. Also, similar to the 
studies done in references 11 and 13 for the midsection of elbow, this study was further extended to 
understand the stiffening effect on stresses due to pressure at the tangent point. The main objective 
of this study was to identify the importance of selecting proper locations for mounting strain gauges 
and utilizing the obtained results to anchor dynamic models for accurate structural and life 
assessments of the SSME ducts under dynamic environment. The finite element method was utilized 
in this study as used in references 9, 11, and 12. 


II. FINITE ELEMENT MODEL 


The finite element models of a 90° elbow with ends terminated by tangent pipes were con- 
structed using the ANSYS isoparametric quadrilateral shell elements. 10 The geometric dimensions, 
boundary conditions, and loadings were taken from a model analyzed in the literature 11-13 to assure 
the validity of a finite element model used in this study. The basic assumptions of homogeneous, 
isotropic thin shells within the range of linear elastic behavior were also used in this study. The 
models have a ratio of R/r , r/t, and pipe factor ( Rt/r 2 ) of 3.0, 29.13, and 0.103, respectively, where R, 
r, and t represent the bend radius of curved pipe, mean cross-section radius of curved pipe, and 



pipe- wall thickness, respectively. In the model, an element aspect ratio was approximately 2:1 for 
tangent pipes and 1:1 for elbow. Figure 1 illustrates the geometry, loading axes, and sign 
conventions for the analysis models. A three-dimensional finite element model contains 1,408 shell 
elements and 4,208 nodes with six degrees of freedom at each node shown in figures 2 and 3. One 
end of the tangent was fixed, and nodal forces for each loading case were applied at the other end. 
The moment (1.737xl0 6 in-lb) used in references 11 to 13 was taken. In addition, internal hoop 
stress pressure (p — 26.66 lb/in^) 11 equivalent to the net pressure (800 Ib/in^) 13 obtained from the 
relationship 15.0xp/0.5 = 800 lb/in^) was also applied to all elements to simulate the pressure 
stiffening effect. For simplicity, nominal bending moments were applied by means of a statically 
equivalent, linearly varying, axial nodal loading pattern as used in references 14 and 15. Similarly, 
torsional moment was applied by means of a statically equivalent tangential nodal loading 
pattern. 14 15 In the area of the tangent pipe, away from the applied loads, nominal bending and 
torsional stresses calculated by the fundamental formula, Mr/I for bending and Mr/J for torsion, were 
observed, where 7 and 7 are a second moment of inertia and a polar moment of inertia, of pipe cross 
section, respectively. Also, it was observed that the stress profiles due to in-plane bending moment 
(figs. 4 and 5) at elbow midsection were similar to those given in references 11 and 13. These 
observations verified the validity of the finite element model and loading application technique used 
in this study. 


in. EFFECT OF TYPES OF LOADS 


The nodal stresses in the longitudinal and circumferential directions on the outside and inside 
shell surfaces were computed for each loading condition. From these computed stresses, the non- 
dimensional stress intensification factors along the cross sections of the duct were calculated. A 
stress intensification factor is defined as the ratio of calculated stress intensity at a point to the 
nominal stress in a piping component. Nominal stress is defined as Mr/I and Mr/J for bending and 
torsion, respectively. From the analyses, the computed circumferential (hoop) and longitudinal 
(axial) stresses, in terms of the stress intensification factors, at section BB (fig. 1) are presented in 
figures 4 through 10. The results/observations for each loading condition are discussed in the 
following subsections. 

3.1 In-Plane Bending Moment ( M z ) 

The longitudinal and circumferential nodal stresses on the inside and outside surfaces at sec- 
tions AA and BB were computed. As shown in figures 4 and 5 for section AA, the maximum circum- 
ferential and longitudinal stresses occurred on the inside surface at 0° and outside surface at 11.25°, 
respectively. The circumferential stress factor was larger than the longitudinal stress factor. Also, in 
this case the finite element analysis results were compared with those derived in references 11, 13, 
and 16, and good agreements were obtained. Whereas, at section BB shown in figures 6 and 7, the 
maximum circumferential and longitudinal stresses occurred on the inside surfaces at 0° and 22.5 , 
respectively. It was also observed that at this section the maximum longitudinal stress was larger 
than the maximum circumferential stress. 3 As shown in table 1, the position angles of maximum cir- 
cumferential stresses were the same between sections AA and BB. However, the position angles of 
the maximum longitudinal stresses were different between sections AA and BB. These observations 
highlight the importance of selecting the proper locations for the mounting of the strain gauges at the 
tangent point. 
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3.2 Out-of-Plane Bending Moment ( M y ) 

At section BB, as shown in figures 8 and 9, it was found that the maximum circumferential 
and longitudinal stresses occurred on the outside surface at -45° and on the inside surface at 
-11.25°, respectively. Similarly to in-plane bending moment, the maximum longitudinal stress factor 
was larger than the maximum circumferential stress factor. However, comparison of figures 6 and 7 
with figures 8 and 9 identified that the maximum circumferential stress factor for out-of-plane 
bending moment was larger than that of in-plane bending moment. However, the maximum 
longitudinal stress factor of out-of-plane bending moment was smaller than that of in-plane bending 
moment. This comparison shows that the results obtained at section BB are contrary to those 
obtained in section A A as given in reference 16. 

3.3 Torque ( M x ) 

Similarly to the out-of-plane bending moment, since a pure torsional moment exits only on 
one tangent pipe at the end of the elbow, any additional steps are not needed to calculate the stress 
factors at section BB. In this load case, the maximum shear stress intensification factor was 
observed on the outside surface at 0° as shown in figure 10. In general, the stress factor obtained 
from a finite element model agreed well with the currently used value given in the ASME section in 
code. 17 

3.4 Effect of Internal Pressure (M z +p) 

The stiffening effect of internal pressure on the flexibility and stress factors of elbow mid- 
section has been studied and reported in references 4, 6, 11, 12, and 13. However, as mentioned 
earlier, these studies did not address the effect of bend stiffening at the tangent point. In this study, 
the resulting stresses were compared to those of cases subjected to moment loads alone to see the 
effect of internal pressure. However, only the in-plane bending moment case is described since the 
conclusive results were not observed from the out-of-plane bending and torsional moment cases. At 
section AA, the finite element analysis results showed an extremely good agreement with the test 
results given in references 11 and 13. Although this combined load produced lower stress factors 
compared to those of in-plane bending moment alone, the basic stress profiles were similar at 
sections AA and BB for both the cases as shown in figures 4 through 7. From this observation, it 
was concluded that internal pressure stiffening effect at section BB also results in reduced bending 
stress at section AA. 13 The position angle of maximum stress components is also shown in table 1. 
As shown in this table, contrary to section AA, the maximum longitudinal stress at section BB 
occurred on the outside surface. 


IV. CONCLUSIONS 


In general, a design of duct systems does not require a detailed knowledge of stress profile 
along the cross section at the tangent point. However, for the duct system in space vehicles a clear 
understanding of stress profiles at the weld location (tangent point) is required since it operates 
under the harsh environment and large dynamic loads. Additionally, it does not have the liberty to be 
over-designed because of the weight penalty. Therefore, this study was conducted using the finite 
element models and, consequently, the detailed stress profiles along the cross section at the tangent 
point of a 90° elbow subjected to various types of loads were obtained. 
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The results show that a lack of the knowledge of detailed stress profile and peak stress 
location (position angle) for the different types of loads might result in erroneous structural and life 
predictions. In general, since the dynamic models used to predict the loads are anchored with the 
strain gauges data obtained from the “hot-fired” systems, the observations of the position angles of 
maximum stress components are believed to be the most useful finding from this study. In other 
words, it means that the strain gauges mounted, in general, at 0, 90, 180, and 270° locations will 
underpredict the strains which may overpredict the fatigue life. Also, it is worth mentioning here that 
although the 30-in diameter duct was modeled in this study, the results obtained here are basically 
applicable for other 90° elbows terminated by tangent pipes with similar ratio(s) of R/r, r/t, and 
Rt/r 2 . 
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Table 1. Position angles of maximum stress components at sections AA and BB. 
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Figure 1. Geometric dimensions of pipe elbow and loading axes 
(a, b, c, d, e, and f indicate the positions for cross section). 
































Figure 2. Front view of Finite element model. 
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Figure 4. Circumferential stress factors of in-plane bending moment 
with (w)/without (wt) internal pressure at section AA. 



Figure 5. Longitudinal stress factors of in-plane bending moment 
with (w)/without (wt) internal pressure at section AA. 
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stress intensification factors 



Figure 6. Circumferential stress factors of in-plane bending moment 
with (w)/with (wt) internal pressure at section BB. 



Figure 7. Longitudinal stress factors of in-plane bending moment 
with (w)/without (wt) internal pressure at section BB. 




stress intensification factors 



Figure 8. Circumferential stress factors of out-of-plane bending moment 
without (wt) internal pressure at section BB. 



Figure 9. Longitudinal stress factors of out-of-plane bending moment 
without (wt) internal pressure at section BB. 
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